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Abstract 

We study the formation of stationary localized states using the discrete non- 
linear Schrodinger equation in a Cayley tree with connectivity K. Two cases, 
namely, a dimeric power law nonlinear impurity and a fully nonlinear sys- 
tem are considered. We introduce a transformation which reduces the Cayley 
tree into an one dimensional chain with a bond defect. The hopping matrix 
element between the impurity sites is reduced by Xj^jK. The transformed 
system is also shown to yield tight binding Green's function of the Cayley 
tree. The dimeric ansatz is used to find the reduced Hamiltonian of the sys- 
tem. Stationary localized states are found from the fixed point equations of 
the Hamiltonian of the reduced dynamical system. We discuss the existence 
of different kinds of localized states. We have also analyzed the formation 
of localized states in one dimensional system with a bond defect and nonlin- 
earity which does not correspond to a Cayley tree. Stability of the states is 
discussed and stability diagram is presented for few cases. In all cases the 

total phase diagram for localized states have been presented. 
PACS numbers : 71.55.-i, 72.10.Fk 
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I. INTRODUCTION 



One well studied nonlinear equation in condensed matter physics and optics is the dis- 
crete nonlinear Schrodinger equation (DNLSE). The DNLSE is a nonintegrable standard 
discretization of the integrable nonlinear Schrodinger equation [|^. The DNLSE in one di- 
mension in its general form is a set of n coupled nonlinear differential equations. 

dCm 

where Kn,m+i = K^+i,™; and m = 1, 2, 3, ...n. (1) 

In Eq.(l) the nonlinearity appears through functions fm{\ Cm |) and Xm is the nonlinearity 
parameter associated with the m-th grid point. Since, J2m I Cm P is made unity by choosing 
appropriate initial conditions, | Cm P can be interpreted as the probability of finding a 
particle at the m-th grid point. The analytical solutions of Eq.(l) in general are not known. 
Numerous works, both analytical as well as numerical, on the DNLSE have been reported 
PI^T^. As for its application particularly in condensed matter physics, we cite among oth- 



ers the exciton propagation in Holstein molecular crystal chain 0. In general, the exciton 
propagation in quasi one dimensional systems |]T5| having short range electron phonon inter- 



action can be adequately modeled by the DNLSE. Other examples include nonlinear optical 
responses in superlattices formed by dielectric or magnetic slabs |]TB| and the mean field 
theory of a periodic array of twinning planes in the high Tc superconductors |]T7 |. 



One important feature of the DNLSE is that this can yield stationary localized ( SL ) 
states. These SL states might play a relevant role in the nonlinear DNA dynamics [0 and 



also in the energy localization in nonlinear lattices |jT9|]. It has been shown that the presence 



of a nonlinear impurity can produce SL states in one, two and three dimension p0|^8|. The 



formation of SL states due to the presence of a single and a dimeric nonlinear impurity in few 



linear hosts has been studied in details |24]. The same problem is also studied starting from 



an appropriate Hamiltonian |^7[ . The fixed point of the Hamiltonian p5|-[28[] which generates 



the appropriate DNLSE can also produce the correct equations governing the formation of 
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SL states. We further note that the appropriate ansatz for the dimer problem has been 
obtained in our earher analysis |2^. Furthermore, the formation of inter site peaked and 
dipped stationary localized states has been studied using the dimeric ansatz. The effect of 
one nonlinear impurity as well as a dimeric impurity in an otherwise perfect nonlinear chain 
on the formation of SL states has also been studied. The other important aspect to be 
considered is the effect of connectivity on the formation of SL states. The simplest way to 
study this is to consider the Cayley tree. The effect of one power law nonlinear impurity in 
an otherwise perfect Cayley tree has already been considered. |Q So, we plan here to study 
the effect of a dimeric nonlinear impurity in an otherwise perfect Cayley tree. We further 
consider the fully nonlinear Cayley tree. For this purpose a transformation is devised to map 
the system to an effective one dimensional chain with a bond defect. Consequently, from 
the symmetry consideration, the dimeric ansatz is found to be suitable for this study. The 
appropriate form of the monomeric ansatz is also found for this case to study the formation 
of on-site peaked solitons. 

The organization of the paper is as follows. In sec. II we introduce a transformation 
which reduces the Cayley tree into a one dimensional system with a bond defect. This 
transformation has also been checked through Green's function analysis in Appendix I. SL 
states in the Cayley tree in presence of a dimeric nonlinear impurity is discussed in sec. III. 
In sec. IV, fully nonlinear Cayley tree is discussed. Finally in Sec. V we summarize our 
investigations. 



II. TRANSFORMATION OF CAYLEY TREE TO ONE DIMENSIONAL SYSTEM 

The structure of a Cayley tree with connectivity, K = Z ~ 1 is shown in Fig.(l). Z is 
the coordination number. We pick up a connection and its two ends are numbered and 1 
respectively without any loss of generality. Furthermore, all points in a given generation lie 
in a shell. Shells are number by n and n & Z as shown in Fig.(l). In a perfect Cayley tree 
the number of points in the n*^ shell is K^~^ if n > 1 and i^T'"' if n < 0. We further note 
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that for a perfect Cayley tree all points in a given shell have identical neighborhood. 

We consider now the motion of a particle on a Cayley tree with connectivity, K. In the 
tight binding formalism with nearest neighbor hopping only equations governing the motion 
of the particle are 

dC ^ ^ _ 

— = KCn-\-i + Cn-i + e„C„, n > 1 



dt 

. dCn 

i—r- 
dt 

dCi 
I—— 

dt 

.dCo 



= KC_\n\-l + C'-|n|+l + CnCn, n < 

= KC2 + Co + tiCi 

= + Ci + eoCo (2) 



dt 

Here denotes the probability amplitude at any point in the n*'* shell and all points in the 
n^^ cell have the same probabihty amphtude because of identical neighborhood. The nearest 
neighbor hopping matrix, V has been taken to be unity without any loss of generality. It is 
also assumed that all points in a given shell arising due to a specific organization have same 
site energy. We, however, note that in our subsequent work with the DNLSE this assumption 

will be automatically satisfied. The normalization condition for the site amplitudes gives 

00 

K 



YlK\^^\\Cn\' + ^Y.K''\Cn? = l. (3) 



n=l 

We now make the following transformations, (i) r = VKt, (ii) e„ = en/VK, (iii) Cn = 
K-i^-^y^Cn, forn > 1 and (iv) = i^-l^l/^C^ for n < 0. After substituting these 

transformations in Eq.(2) we finally obtain 
.dCn 



dr 



= Cn+i + Cn-i + enCn, for n> 1 and n < 0. 



— — C2 -\ — t^Cq + eiCi, 
dr JK 



•"^^^ C_i + 4^Ci + eoCo. (4) 



dr 4k 

Furthermore from Eq.(3) normalization condition reduces to Y^^oo iC'nl^ = 1- So, the motion 
of a particle on a Cayley tree is mapped to that on a one dimensional chain. However, in 
this chain the nearest neighbor hopping matrix element between the zeroth and first site 



is reduced from unity to I/VK. In appendix I we show that the Green's function Gofl{E) 
calculated from Eq.(4) will yield the Gofi{E = E\/K) of a Cayley tree with connectivity K. 

Since we are interested in the DNLSE with general power law nonlinear impurity, in 
our case e„ = e„/v^ = Xn^~^""^^"/^^~^/^|C„|" for n > 1 and e_|„| = e_|„|//^ = 
X-|n|-^ '"'"^^^-^ "'^^^IC'-lnlT > 0. We further note that a model derivation of power 



law nonlinearity is given in ref. When all points have the same nonlinearity strength, 
we have Xn = Xi ^ ^ ^ ■ the other hand for a dimeric nonlinear impurity, Xn = 
xi^nfi + Sn,i),n G Z. Furthermore, the classical Hamiltonian which can generate Eq.(4) is 

oo 1 oo 1 

^= E ^|Cnr-'' + 2 ^ {CnC:^i + C:C^+i) + -VoiCoCt + C,C'o) (5) 

n=— oo ' ^ ^ n=— oo 

where Vq = - 1). 



III. A DIMERIC NONLINEAR IMPURITY IN THE CAYLEY TREE 

Here we are interested in the possible solutions for SL states due to a dimeric impurity. 
Hence we assume that 

where 

0m = {sgn.{E)r])"^^^ 01 for m > 1 
and 

0_H = (sgn(^)77)l™l0o form<0. (6) 

Eq. (6) is the exact form of 0m in the presence of a dimeric impurity and can be derived from 
Green's function analysis, rj e [0, 1] is given by r/ = (jE\ — \^E'^ ~ Aj /2. Sgn(ii^) denotes 
the signature of E. We further define jS = 0o/0i if |0i| < |0o|- Otherwise we invert the 
definition of (3. Because of the symmetry in the system we will obtain the same result. So, 
apparently /? G [—1, 1]. However, we shall show later that for x > 0, negative values of P 
except (3 = —1 are not permissible. The analytical argument showing the impossibility of 
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such a situation is presented in ref. p4[. Of course, negative values of /3 > — 1 will produce 
SL states in the antisymmetric set if x < 0. So, for x > 0, /3 G [0, 1] and for x < 0, 
(3 G [—1,0]. This assertion will also be substantiated here in due course. Now, from the 
normalization condition we get 



Introducing Eqs.(6) and (7), and the definition of j3 in the Hamiltonian (Eq.(5)) we get an 
effective Hamiltonian, H^ff where 

i/e// - sgn(E)r; + ^ ^ + ^ Jl^WW^' 

and V = for the Cayley tree otherwise it is a parameter. The Hamiltonian consists of 
two variables, namely, (3 and r] and two constants, namely, x ^i^nd a. The stationary localized 
states correspond to fixed points of the reduced dynamical system described by H^jf. 

1. |/3| = 1 

We here consider the case where \(3\ = 1. For x > 0, /3 = 1 yields the symmetric set, 
while for /3 = — 1 we get the antisymmetric set. In this limit the relevant equation governing 
the formation of SL states is obtained by setting dHeff/drj = 0. This in turn yields 

X sgniE) - sgnif3)Vv 

Case 1: a = 

Here we have a linear dimeric impurity in a Cayley tree. Since we are considering x > 0; fo^^ 
the symmetric case we must have sgn(£') = 1. On the other hand for antisymmetric case 
sgn(£') can take both the signs. If \^ < 1, for the symmetric case F{0,r]) has a divergence 
at ?7 = > 1. Since r] e [0, 1], this divergence occurs beyond the permissible range of rj. 
However, in the permissible range of r/, ^ = (1 — Vi])^'^ > 0. So, -F(0, i]) is a monotonically 
increasing function of rj and it assumes the permissible maximum value at ?7 = 1. This in 
turn then yields Xcr = (1 ^ V") and a SL state will be obtained if x ^ Xcr- For the dimer in a 



Cayley tree we then need x ^ Xcr = ^/K — 1. On the other hand for > 1, the divergence 
at ?7 = ^ is in the permissible range of 77. So, even with an infinitesimally small value of x, 
we shall obtain a SL state. However, for x > 0, 77 G [0, 

We consider now the antisymmetric case with sgn(i?) = +1. It is easy to see that -F(0, 77) 
is a monotonically increasing function of rj. So, -F(0,r7) takes the maximum possible value 
at 7] = 1. This in turn gives Xcr = 1 + 1^. ^ is implicitly assumed to be positive. So, we 
shall get a SL state for the Cayley tree if x ^ Xcr = \/K + 1- For sgn(i?) = —1, F{0,ri) 
will diverge at 77 = ^. So, if < 1, we shall not obtain any SL state below the band of 
the host system. On the other hand if V > 1, F{0,ri) diverges in the permissible range of 
rj. However, we also have F(0, 1) = Thus, Xcr = V — 1 and a SL state below the band 
will be obtained if x ^ Xcr = V — 1. We now summarize our findings on the linear dimer. 
(i) V < 1. No SL state will be obtained if x < 1 ~ V". There is one SL state for {1 — V) < 
X< (l + V). But there are two SL states if X > (l + V) p9|. (ii) ^ > 1. IfO<x< (^-1), 
there are two SL states. One appears above the band and the other lies below the band. If 
{V — 1) < X < (1 + W6 have one SL state above the band. For x > (l^ + 1), we get two 
SL states and both appear above the band. 
Case II: a 0. 

Here we consider two cases, namely, < 1 and V > 1 separately. 
(A) V < 1. Since x is taken to be positive, in Eq.(9) we need sgn{E) = sgn(/9) = +1. Again 
the divergence of -F(cr, 77) at 77 = is of no consequence. Furthermore, -^(0", rj) has at least 
one maximum at 1]^ G [0, 1]. In fact, F{a,(3) has only one maximum. So, there will be a 
Xcr so that for x > Xcr ^6 shall obtain two SL states and for x < Xcr; there will be no SL 
state. On the other hand in the antisymmetric case we have sgn(/3) = —1 but sgn(£') can be 
either +1 or -1. In the first case (sgn(_E')=+l), F{a,ri) has no divergence but F(cr, 0) = = 
F(cr, 1). So, F{a,ri) has at least one (actually one) maximum at r]'^ e [0,1]. Consequently, 
we shall get another critical value of x, say Xcr so that if x > Xcr shall obtain two SL 
states. For x < Xcr there will be no SL state. We further note that Xcr > Xcr- the second 
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case (sgn(£') = -1), F{a,rf) diverges at 77 = ^ > 1. Furthermore, F{a,r]) should be positive. 
Hence, r] > y. Since allowed values of 77 G [0, 1], no SL state will be obtained below the 
host band. So, we shall get three regions having no, two and four SL states. Furthermore, 
in half of the states 77 ^ 1 as x — > 00. So, these are unstable states. Equations for critical 
lines in the (x, c) plane separating three regions are given in Appendix II. 
(B) > 1. In the symmetric case F{a, rj) diverges at 77 = < 1 and F((T, 0) = = F(cr, 1). 
Furthermore, F{a,r]) > for 77 < ^ and F{a,r]) < for 77 > ^. Hence, a SL state will 
always be obtained for x > 0. The maximum value of rj the SL state can take is ^ and 
this will happen if x 0. In the antisymmetric case when sgn(£') = +1, F{a,r]) has no 
divergence for 77 e [0, 1]. Since, F{a, 0) = = F{a, 1), F{a, rj) has a maximum at 77^ e [0, 1]. 
So, there will be a critical value of x such that x < Xcr no SL state will be obtained. On 
the other hand for x > Xcr, we shall have two SL states and in one of these states 7; — > 1 as 
X — > 00. So, one of the states is an unstable SL state. Equation for the critical hne in (x, cr) 
plane separating these two regions is also given in Appendix II. For sgn{E) — —1, F{a,r]) 
again diverges at 7; = ^. However, for F{a,rj) to be positive we need ^ < < 1. Since 
F((7, 1) = 0, we shall always get a SL state irrespective of the value of x- Furthermore, the 
SL state will appear below the host band. The minimum value of rj that a SL state in this 
case can attain is y. This will happen if x ~ and rj ^ 1 as x — 00. So, this is an unstable 
SL state. Finally for y > 1, we have two regions. The region below the critical hne has two 
SL states and that above the critical hne has four SL states. As usual, half of these states 
are unstable. 

2. \(3\^1 

We consider here the scenario, \P\ ^ 1. -ffe// (Eq.(8)) has now two dynamical variables, 77 
and (3. So, the relevant equation governing the formation of SL states is obtained by setting 

dHeff /dXi = 0, where Xi = r] and X2 = [3. We note in passing that for a = 0, dHgff/d(3 
— yields P — ±1. From the first condition {dHeff/dr) — 0) we obtain 
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So, 77 is a symmetric function of (5 and (5~^ as enunciated earlier. Furthermore, since 77 > 0, 
if y > 0, sgn(£')=sgn(/9). In subsequent discussion we assume that sgn(/3) =+1. When 
\j3\ — >• 1, from Eq.(lO) we obtain, 

Vu^-^— (11) 

The second condition {dHgff/d(3 = 0) yields 

1 _ sgn(/3)|/?|(l-|/3r) . ,y/2 

X y(l + /32)'^/2(l-/32) V I J ^ ^ 

Since (1 — 77^) for 77 e [0, 1] is a positive semidefinite quantity, for > 0, x and (3 will 
possess the same sign. In other words, SL states with ^ 1 in the antisymmetric set 
are not possible. It is trivially seen that the right hand side of Eq.(12) is also a symmetric 
function of /3 and We further note that introduction of Eq.(lO) in Eq.(12) makes the 
right hand side of Eq.(12) an explicit function of a and p. We call this function F{a, (3). 

Since cr = 2 is physically more relevant we consider this case in detail. In this situation 
we have 

- = r,{l-r,^)=F{2,r,)=g{r,) (13) 
X 

We note that g{ri) has one and only one maximum at 7;^ = | for r] € [0, 1]. Furthermore, 
riu — {2V)~^. So, if < 0.5 or X > 4, 77^ > 1. Since 51(77) in this case is defined for 
77 e [0, 1], the maximum of g{r]) lies in the permissible range of 77. On the other hand if 

V > 0.5 {K < 4), 77„ < 1. So, rj G [0,77j. Then for rjm to stay in the allowed range of r], 
we need 77^ < 77^. This in turn yields V < 0.8660 or K > 1.33. So, there will be a lower 
critical value of x, Xcri such that x < Xcri there will be no SL state and for x > Xcri there 
can be two SL states. From Eq.(13) we further obtain Xcri = Xch/VK = 2.5980. Again for 

V < 0.5, since 77 e [0, 1] and g{0) = = ^(1), we shall obtain two SL states for V < 0.5 and 
X > 2.5980. In one of the states 7/ ^ 1 as x ^ oo- So, one state is an unstable state. On 
the other hand, 0.5 < V < 0.8660, we get Xcm = if^zi and x > Xcm, we get only one SL 
state. Then, if < 0.8660 and Xcri < X ^ Xcm we have two SL states. In one state 77 ^ 77^ 



as X — ^ oo. So, one state is unstable. Thirdly, for V > 0.8660, r] G [0,r]u] and r]u < r/m- So, 
g{r]) takes the maximum value at rj^ and the corresponding critical value of x is Xcm- We 
note that for = 1, Xcm = 8/3 p^j27[] and for V = \f2, Xcm = 3.232. For x > Xcm we shall 



obtain one stable SL state. 

We now consider the general a. Substituting = 1 in Eq.(ll) we obtain a' = 2V/{1 — V) 
if V < 1. For the Cayley tree it translates to cr' = 2/(v^ ~ simplicity we break the 

discussion in two cases. 
Case I: a > a' 

Here rju > 1. So, r] G [0,1] but j3 G [0,/3n] where Pu < 1- Furthermore, F{Pu,0') = = 
F(0, cr). So, for a given a > a', F{P,a) has at least one maximum at Pm ^ [O,/^^]. For 
K = A or V = 0.5, a' = 2. In Fig.(2) we have plotted F{/3,a) for a = 2.5. We see that 
there is only one maximum. Due to the maximum at /5m, in the (x, a) plane there will be 
a critical line separating the no state region from the region containing two SL states. The 
equation of the critical line is xir^ = [F{Pm, cr)]^^ where Prn is the solution of dF/dP=0. 
Since in one of the states 77 ^ 1 as x ^ 00, the two states region has one unstable state. 
Furthermore, as we go along the critical line /? — 0, 77 and x ^ oo- This implies that 
in the stable state, the amplitude gets preferentially localized in one of the dimer sites as 
cr —i> 00. X ^ 00 because in the limit a monomer localized state is formed. 
Case II: a < a' 

Since ?7„ < 1, /3 G [0, 1]. Consequently, there will be a critical value of x, Xc? given by 

For X > Xc?) W6 shall obtain at least one SL state. Since Xc? ^ 00 as cr ^ and a ^ a' 
{r]u 1), Xc? "wiii assume a minimum value at amin- crmin is obtained from dx^cr /da =0. 
However, it gives a very complicated algebraic equation in a. When a — > 0, 0. So, 
we obtain a SL state localized mostly on the dimer. Precisely for this x^cr ^ C)0 as a ^ 0. 
Again a increases, rj increases. This will require lower values of x- So, for < a < (Jmin and 
X ^ xfr 1 the system yields a SL state with /5 7^ 1. On the other hand, for amin < o" < cr', 

10 



albeit i] increases towards unity as o" ^ cr', Xc? increases towards infinity. So, the SL state 
obtained for a > Omin (in fact cr„. defined later) and x ^ X^cr unstable. 

F{j3,a) also develops a local maximum at jSm < 1- This is shown in Fig. (2) for = 4 
(1^=0.5) and a = 1.5. So, there will be a critical value of a, Ocr defined as follows. If 
(J = (Jcr + 5 and 5^0, there exists an e — > depending on 5 so that /?m = 1 — e- So, 
for a > there will be a lower critical value of Xc?'' and x < xir"* state will 

be obtained. It is trivially seen that xir"* — iPm, <^)]~^ ■ For a given a > the upper 
critical value is Xc?(^)- Then in the (x, o") plane we have a V-shaped region with boundaries 
acr < cr < cr' and x > xi^H'^cr) but Xcr'*!^) > X > xir''!^)- This region contains two SL 
states with /5 7^ 1. However one of the states is unstable because 77 — r^^ as x 00. We 
further note that the lower boundary of the V-shaped region joins smoothly with the critical 
line for a > a' . 
Case III: V > 1 



The scenario is very similar to the V = 1 case discussed in ref. [Q. However acr moves 
to a higher value depending on the magnitude of V. 

The total phase diagram of SL states for V = 0.5 {K = 4) is shown in Fig. (3). The 
number of possible SL states in each region is indicated. We again note that the maximum 
number of SL states is six. The total phase diagram for V = \p2 is shown in Fig. (4). The 
maximum number of SL states is found to be five. Furthermore, the a = 2 line has no 
special significance. The stability diagram of SL states for dimeric nonlinear impurity 
in a Cayley tree is shown in the (77, a) plane in Fig. (5). 



IV. ALL SITES NONLINEAR 

1. Inter site peaked and dipped solutions 

We consider here the formation of SL states in a Cayley tree with each site having a power 
law nonlinearity and the same coupling constant, x- But due to the transformation proposed 
here, the problem reduces to the study of formation of SL states in a one dimensional system 
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with a bond defect between the zeroth and the first site. The tunnehng matrix between the 
sites, V is reduced by a factor of l/y/K and site energies are 

en = V'-WCnV if n > 1 
and 

6_|„| = V^I"lx|C!|^„| if n > 0. (15) 

Since we have a bond defect, from the symmetry consideration, we use the dimer ansatz 
first to study the problem. The apphcation of the method outhned here and in other places 
^,27] yields the equation governing the formation of SL states in this system. It is given 



by 

2-/2 ^(1 _ rfyn _ (y^Y^ 



X (sgn(E) - sgn(/3)V^r/) (1 - r^'^iVr^Y) 



G(n,a) (16) 



rj is defined earlier [^. We are also assuming that x > 0. Since [1 — (VriY] —a IniVr]) 
as cr ^ 0, G{rj, a) ^ as a ^ 0. Consequently x ^ oo. This implies that no SL state will 
be formed in this limit. We now consider various cases. 

(I) Symmetric Case 

Here sgn(£') = sgn(/5) = +1. We note that G{r],cr) have a removable singularity and 
a divergence at rjo = 1/V and rji = l/V^ respectively. But for V < 1, rji > 1. So, 
the singularity of G{r], a) at r/i will not play any role in the formation of SL states. Again 
(7(0,0") = = G{l,cr). So, G{ri,a) will have at least one maximum at 1]^ G [0,1]- It can 
be seen numerically that G{r],a) has only one maximum for rj G [0,1]. Consequently in 
the {x, 0") plane there will be a critical line separating two states region from the no state 
region. Since one of the states in the two states region — 1 as x ^ c>o, it is an unstable 
state. On the other hand for V > 1, G{ri,a) has a divergence at r/i < 1. So, in this case 
the system will always produce a SL state even if x is infinitesimally small. Furthermore, 
lim^^oG{rii — e,a) ^ oo from the positive side only. So, for x > 0, rjmax = — > as 
X ^ oo. Hence this is a stable SL state. 

(II) Antisymmetric Case 
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In this limit sgn(/3) = -1 but sgn{E) can be either +1 or -1. We note that for sgn(i?)=+l, 
G{rjo,a)=0 here. But for ^ < 1 both rjQ and rji lie beyond unity. Since G{ri,a)=0 both at 
ri=0 and //=!, it has a maximum at rim G [0, 1]. This imphes that in the (x, cr) plane there 
will be a critical line. This line again separates the no state region and the two states region. 
Furthermore, in the two state region one state will be unstable for the same argument given 
earlier. 

For > 1, we note that G{r],a) goes to zero and infinity at t]q and r]i respectively. 
Furthermore, if cr is finite, rjo < rji. So, for r] G (?7o,r7i), G{rj,a) is negative. This in 
turn implies that lim^^QG{rji — e, a) — > — oo and lim^^QGijji + e, a) — > oo. Consequently, 
lim^^QG{l — e, 0") — s> from the positive direction for 77 G [?7i, !]• Therefore, we shall always 
obtain a SL state even if x is infinitesimally small. In this SL state, however, rjmin = Vi 
as ^ 1, X — i> 00. So, this is an unstable state. We further note that G(0, a)=0=G{riQ, a). 
Then there will be a maximum of G{r],a) at rim G [0,?7o]. So, there will be a critical line 
in the (x, cr) plane also. This line will separate one state region and three states region. 
It is further seen that two of the states in the later region are unstable. We can also have 
sgn(£')=-l. However, for < 1 no SL states will be obtained in this limit. On the other 
hand for V > 1, lim^^QG{rji =F c") ^00. So, we shall always get a SL state below the 
band even if x is infinitesimally small. However, this SL state is unstable. 

We now combine our results to obtain the phase diagram. For < 1, we have three 
regions, namely I, II and HI containing no SL state, two SL states and four SL states 
respectively. This is shown in Fig. (6) for K = 4 or V = 0.5. The stability diagram of 



SL states for this case is shown in Fig. (7). However, for V > 1 we do not have no SL state 
region. Instead we have a three state region separated from a five state region by a critical 
line. One of these states appears below the band. In the three state region we have two 
unstable states while in the other region we have three unstable states. For V = a/2, the 
phase diagram is shown in Fig. (8). 
2. On site peaked soliton 

We discuss here the formation of on-site soliton in the transformed system described by 
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the system given by the Eq.(2). For this purpose we first consider a power law nonhnear 
impurity with strength, x embedded at the zeroth site. After introducing the dimeric ansatz 
in the appropriate form of the Hamihonian we obtain 

- 7T-2 (ttI)"'"" + =^"<^' " + (rrl) (1^' 

Again relevant equations are obtained by dHf.ff/dXi=0 where Xi = f3 and X2 = rj. After 
a trite algebra we then obtain (3 = sgn(E)V"?7 and 

= ^(^-^)^^' = /(, a) (18) 

X (i + yv)"/'(i-^V) 

For V = Eq.(18) describes the formation of SL states due to a nonlinear impurity 



in a Cayley tree. This has been discussed in detail in ref. |2J]. So, we see that the dimeric 
ansatz reduces to the appropriate monomeric ansatz. When > 1, f{ri, a) has a divergence 
at ?7„ = l/V. Furthermore we have, /(0,(t) = = /(l,cr) and lim^-^o /(^m — o") oo- 
So, we shall obtain two SL states even if x is infinitesimally small and cr > 0. However, one 
state will appear below the band. In this state rjmin = Vu and rj ^ 1 as x ^ 00. So, this is 
an unstable state. For, cr = 0, this state will appear if < x < (^^ — !)■ 

To study the formation of on-site peaked SL states in the fully nonlinear chain we put 
f3 = sgn{E)VT] in H^ff given by Eq.(17). We then set dHeff/drj = 0. After a trite algebra 
we finally obtain 

1 r]{l-r]^y^^ {1 + V''+^7]''+^){1-V''t]'') 



X (1 + V^ri^y^^ (1 - V-ri-+^ f (1 - V^T]^) 



fi{v,cT) (19) 



When V=l, Eq.(19) reduces to the relevant equation in ref. |2^. We note that fi{ri,a) has 

— a 

a removable singularity at rju and a divergence at r^i = V <^+2 . However, for V < 1 the diver- 
gence at rji does not play any role in the formation of SL states. Since /i(0, cr)=0=/i(l, a) 
we expect at least one maximum of f\{ri,a) at rjm ^ [0, !]• It is seen numerically that for 
rj G [0,1] fi{ri,a) has only one maximum. So in the {x,<^) plane there is a critical line 
separating the no state region from the two states region. Of course, one of the states is 
unstable. In Fig. (6) the critical line is shown by solid curve for K = 4. On the other hand, 
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for y > 1, fi{r], a) diverges at rji. Since, lime-^ofi{r)i ^ e, o") — > oo we have also two states 
and these states are formed even if x is infinitesimally small. One of the states is unstable. 
However, numerical calculation shows that there exists a critical value of a say acr such 
that for (7 > acr there will be a four state region bounded by two critical values of x- For 
example, for V — ^/2 ~ 3.85. But the four states region occurs at larger value of x- Two 
of the states are again unstable. 

V. SUMMARY 

The DNLSE with general power law nonlinearity is used to study the formation of sta- 
tionary localized states in the Cayley tree. Two cases, namely, a dimeric nonlinear impurity 
and the fully nonlinear system are considered. To facilitate the study a transformation is 
devised to map the system to an one dimensional system with a bond defect. We also note 
in passing that the problem can also be mapped to an half infinite chain with a bond defect 
between the zeroth site and the subsequent site. The Cayley tree Green's function for the 
problem can be obtained from the transformed system. This is discussed in Appendix I. 

The formation of SL states is studied by analyzing the fixed point equations of the reduced 
dynamical system. This is obtained by introducing the dimeric ansatz in the appropriate 
Hamiltonian. For the linear dimer our results agree with known results. In case of a nonlinear 
dimer impurity, the system is found to sustain two types of SL states and altogether a 
maximum of six types of SL states can be obtained. In one case the absolute amplitude at 
two sites are equal. In the second category we find states with unequal amplitudes. In this 
aspect our results are very similar to what we obtain for a one dimensional chain. There 
are, however, some differences. In the one dimensional chain the no state region is obtained 
for (T > 2. Furthermore, the V region extends to infinity {a' — > oo). Here, no state region is 
obtained from a = 0, and a' = -^y^Zrj- So, the V region shrinks as K increases. For a = 2, 
we find that for K > 1.33 (V < 0.86), a third critical value of Xcr — Xcr 

jy/K = 2.5980. 

This is a very interesting result. The corresponding Xcr for X = 1 is 8/3. The stability and 
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the phase diagram of SL states are discussed in detail. 

The formation of SL states in the fully nonlinear Cayley tree is also considered. For the 
on site peaked solution, the appropriate ansatz is derived. In the perfect nonlinear chain, a 
three SL states region exists for the on site peaked SL state. For the Cayley tree, this region 
is absent. Instead for all cases, we find a two states region and a no state region separated 
by respective critical lines. Along with this the case where V (in relation to other hopping 
element) > 1 is considered. For this case we show that under certain conditions, states can 
appear both below and above the band. Furthermore, we also find a four states region. 



VI. APPENDIX I 



We show here that a subset of amplitude Green's functions of Eq.(4) with €n — 0, n E Z 
yields the amplitude Green's functions of a particle moving on a Cayley tree with connec- 
tivity, K. We first note that Gn,m{E = E^/K) = -^Cra{^)Gn,m{E). Furthermore, from the 
transformations, we have (5^(0) = Cm{^)/ K^"^-^)/"^ if m > 1 and (5^(0) = C^(0)/is:l'"l/2 if 
m < 0. The Hamiltonian {H) that yields Eq.(4) is H ^ Hq + Hi where, 

Hq^Y^ (a^an+l + Otn+lOtn) 
n 

and 

Hi = (^-^ - 1)^ (4ai + a\ao) (20) 

ttnial) destroys (creates) a particle at the n-th site. We define Vq = {1/\/K — 1). We further 
have G{E) — {E — H)~^ and Gq{E) — {E — Hq)~^. These two operators are related by 
(/ — GqHi)G — Go where / is the identity operator. By defining < n\G{E)\m >— Gn,m{E) 
and with similar definition for Go(n,m) (-E") we obtain from the relation between G and Gq 

Gn,m{E) + Vq (Go{n,l)(E)Go^m(E) + Go(n,0) (-E')G'l,m(-E')) = G'o(n,m) (-E') ■ (21) 

We have then two unknowns, namely, Go,m{E) and Gi^rn{E). After some algebra we obtain 



Go,m{E) ^ 



Gl,m{E) 



1 

D 



1 - K)G'o(i,o)(-E) yoGo(o,o)iE) 
VbGo(o,o) (E) 1 - 14)^0(0,1) (E) J 
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0(0,m) 



(E) 



Go(l,m){E) 



(22) 



and 

D = {l- VoGo^o,i){E)) (l - K)Go(i,o)(i?)) - ^,'^^(1,0) (^) (23) 
We again note that Go(m,n){E) = [sgn{E)f-^+'Goio,o){\E\W~""^ and Go(o,o)(|^|) = 



— 4 for \E\ > 2. Hence, 6*0(0,0) ( I -E"]) = (1 ^ V^)/'']- After some simple algebra 
we obtain D = [{K — l)?7Go(o,o)(|-£'|) + K]/K. Furthermore, from Eq.(22) and relevant 
transformations we obtain 

~ ~ 1 ~ 2K 

Go,o{E) = -j=Go,o{E) = sgn(E) , (24) 

{K-1)\E\ + {K + 1)^Je^-4:K 

Again from Eq.(22) we can easily show that for m > 

^ ) ^^^^ 



|E| + yE'^-AK 

when Go.od-^l) is given by Eq.(24). On the other hand for m < 0, we have 

Go-\m\{E) = -j=^^^j^Go^^\m\{E) = Go^rn{E) (26) 



From Eq.(22) we further obtain, for m > 1 

(\ m-l 
- — =Go,^-m 
\e\ + ^Je^-ak] 

(27) 

On the other hand for m < we obtain 

<5i _|r„|(^) = ^1^1/2 '^i-Ih(-^) = Go,\m\+i{E) (28) 



Since in relation to the Cayley tree we are dealing with a translationally invariant problem, 
the choice of origin is arbitrary. So, (0,1) bond can be rechristened (n, n + 1) without any 
loss of generality. If the shift in the origin is incorporated in the transformation, Go^rn{E) 
and Gi^rn{E) will be transformed to Gn,m{E) and Gn+i,m{E) respectively. Consequently our 
result will agree in full with the calculation in ref. ||29|| . 
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VII. APPENDIX II 



We derive here the equation of the critical hne for the symmetric state with (3=1 and 
for the antisymmetric state. Note that x > and V > 0. The equation to be considered for 
the purpose is 

2-/2 nil-ri^Y^^ 



X sgn(E) - sgn{(3)Vr] 
Then the equation of the critical line in (x, cr) plane is 



F{r],a) (29) 



2^/2 

^ ^) 

In Eq.(30) ± refers to the symmetric and the antisymmetric cases respectively. To find rjf^ 
we set dF/drj = 0. This in turn yields 

aVri^ T ((X + l)r/2 ±1 = (31) 

In Eq.(31) the upper sign refers to the symmetric case and the lower sign to the antisym- 
metric case. For V = 1, we find that 

ari^^ri- 1 = (32) 

Eq. (32) has been derived in ref. [2^ . real positive root of Eq.(31) and it must be 



less than unity. The expression for the r]cr is found to be 



where, Q=^ (A + SaVVW), A=2+6a+6a^+2a^-27a^V^ and B=- 

4-12cr- 12(T^-4(T'^+27(7^V^^. In Eq.(33) Im refers to the imaginary part. While the upper 
sign refers to the symmetric case, lower sign is for the antisymmetric case. 
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FIGURES 



FIG. 1. Cayley tree with connectivity 2. All bonds are of same length. 

FIG. 2. F{(3, a) as a function of /? for /3 7^ 1 in case of a dimeric nonlinear impurity in a Cayley 
tree. Here K=4. Solid, dotted and dashes curves are for cr = 1.5, 2 and 2.5 respectively 

FIG. 3. Total phase diagram of SL states of a Cayley tree in presence of a dimeric nonlinear 
impurity. Here K = 4{V = 0.5). a-nd x" represents the critical lines for symmetric and anti- 
symmetric case respectively. Xcru and Xcrl represents the upper and lower critical line respectively 
for /3 7^ 1 case, amin, <^cr and a' are shown. Numbers indicate the number of possible SL states in 
those regions in the (x, cr) plane 

FIG. 4. Total phase diagram for SL states of a one dimensional chain with a dimeric nonlinear 
impurity and a bond defect in between the impurity sites. Here V = ^2. Numbers indicate the 
number of SL states in those regions. The unnumbered closed triangular small region contains four 
SL states. 

FIG. 5. Stability diagram for SL states in a Cayley tree with a dimeric nonlinear impurity. 
Here K = 4(y = 0.5). Stability of the states in various regions are marked in the figure. 

FIG. 6. Phase diagram for SL states in a fully nonlinear Cayley tree. The solid line is the 
critical line for the on site (zeroth site) peaked solution. The lower dotted line is the critical line 
for the inter site peaked (symmetric) solution and the uppermost line defines the critical line for 
the inter site dipped (antisymmetric) solution. Here K = 4(y = 0.5). For inter site solutions 
region I, II and HI contains no, two and three SL states respectively. For on site solutions there is 
no state below the solid curve and two states above the solid curve. 

FIG. 7. Stability diagram for SL states in a fully nonlinear Cayley tree. Here K = 4{V = 0.5). 
Stability of the states in various regions are marked in the figure. 
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FIG. 8. Total phase diagram for SL states of a fully nonlinear one dimensional chain a dimeric 
nonlinear impurity and a bond defect in between the impurity sites. Here V = ^2. The region I 
contains three SL states and the region II contains five SL states. 
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